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We study a diffusive superconductor-normal metal-superconductor (SNS) junction in an external 
magnetic field. In the limit of a long junction, we find that the properties of such a system depend 
on the width of the junction relative to the length associated with the magnetic field. We compute 
the critical width separating the regime of pure decay (narrow junction) and the regime of damped 
oscillations (wide junction) of the critical current as a function of the magnetic flux through the 
junction. We find an exponential damping of the current, different from the well know Fraunhofcr 
limit which corresponds to the limit of a tunnel junction. In the limit of a wide junction, the 
superconducting pair correlations and the critical current become localized near the border of the 
junction. 

I. INTRODUCTION 

The interplay between magnetism and superconductivity in proximity structures results in a rich physics. The 
experimental observation of the reversal of the critical current (ir coupling) in Josephson junctions with a ferromagnetic 
interlayeri£i£ has recently triggered a lot of interest in the study of the competition between the ferromagnetic and 
the superconducting order j 4 Ferromagnetism has also allowed to generate triplet superconducting pair correlations in 
proximity structures.— Others studies have focused on the Zeeman effect of an external magnetic field in such systems 
(see for example Ref. 0) and on the possibility to obtain a it junction by applying a fields 

In this work, we address the orbital effect of an external magnetic field in a superconductor-normal metal- 
superconductor (SNS) Josephson junction. Most of the experimental applications cited in the previous paragraph 
concern diffusive materials, for which an appropriate formalism has been derived.— In the diffusive regime, the motion 
of the electrons is governed by frequent scattering on impurity atoms: the elastic mean free path l e is much smaller 
than the relevant length scales of the system and the Green function becomes isotropic. In many situations, this 
allows to apply simple one-dimensional models to study the proximity effect. However, including orbital effects in 
the formalism forbids reducing the system to one dimension. As a result, the proximity effect in the presence of the 
orbital effect of the magnetic field has been studied until now essentially numerically or in simple limits (wide and 
short junction or narrow junction) for diffusive hybrid structures*^ and in the clean limits 

It is well established that in the limit of a thin (tunnel) junction the Josephson current changes sign along the 
transverse direction and that the total current exhibits a Fraunhofer-like dependence on the magnetic flux through 
the junction^ Observation of this type of dependence has been extensively used experimentally to confirm the 
Josephson nature of the coupling between superconductors. More recently, shifted Fraunhofer patterns have also 
served as an indicator for the presence of a net magnetization when a ferromagnetic interlayer is usedj 13 ' 14 It has 
been then shown both experimentally^ and theoretically^*^ that in proximity structures, discrepancies from the 
usual Fraunhofer patterns can be present. In particular, the authors of Refs. l9Hl(J have discussed numerically how 
the damped oscillatory behavior (Fraunhofer like) characterizing wide and short junctions is replaced by a monotonic 
exponential decay in narrow junctions. They have also identified the length scale over which the transition between 
the two regimes takes place. 

Motivated by this recent activity and by the rarity of analytical results on the Usadel equation for non one- 
dimensional geometries, we revisit the problem of the diffusive SNS junction in an external magnetic field. We 
consider the limit of a long junction and linearized Usadel equations to obtain analytical results for a two dimensional 
problem. We show that for a narrow junction, the Josephson critical current decays exponentially as a function of 
the flux through the junction. We find the transition point (the critical width of the junction) where this monotonic 
decay is replaced by damped oscillations of the critical current. Finally, in the limit of a wide junction, we find 
damped oscillations with the same period as in the Fraunhofer limit but an exponential decay instead of the purely 
algebraic decay characterizing Fraunhofer patterns. In this regime, the superconducting correlations become localized 
in a small region close to the border of the junction. The method we develop does not rely crucially on the choice of 
particular boundary conditions for the interface between the superconductor and the normal metal: it can be applied 
either to the situation where the SN interfaces are transparent or to systems with finite interface transparency. 

The paper is organized as follows. In Section II we describe the SNS Josephson junction we consider and introduce 
the formalism used throughout the paper. We then compute the superconducting Green function (Section III) and 
the Josephson current (Section IV) for the SNS junction in a transverse field. In Section V we discuss the applicability 
conditions of our method and finally in Section VI we summarize our conclusions. 
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II. SNS JUNCTION IN A TRANSVERSE MAGNETIC FIELD 



We consider a SNS junction in a transverse magnetic field. We introduce a coordinate system with the N layer in 
the yz plane (Fig. 1). The x axis is directed along the junction, and the SN interfaces correspond to the coordinates 
x = 0, L x . The origin of the y axis is chosen in the middle of the N layer and we denote the width L y . The system 
is invariant under translation along the z axis. We take a uniform magnetic field H directed along the z axis and 
neglect the screening of the magnetic field by the Josephson currents^ For simplicity, we further consider that the 
London length is short compared to the length of the junction L x and neglect the penetration of the magnetic field 
in the superconducting electrodes. An exact treatment would require to add twice the London length to L x in order 
to get the total flux through the junction^ 

We assume that the normal layer is strongly disordered and the motion of electrons is diffusive. In this regime, the 
quasiclassical Green functions^ (averaged over the fast Fermi oscillations and the momentum directions) are given 
by the solutions to the Usadel equations^ 

The (nonlinear) Usadel equation in the normal layer takes the form (we follow the conventions used in Ref. [l7l ) 

HDv(gVg) -w[f 3 ,s] =0. (1) 

D = v-pl e /Z is the diffusion constant with l e the elastic mean free path and uj = (2n+ l)TTkgT is the Matsubara 

{ G F \ 

frequency. The Green function g = ( ^ ^ j is a matrix in the Nambu (particle-hole) space, and f a denote the 
Pauli matrices in this space. The gradient operator V contains the vector potential A 

Vs=^_ (v + ^ A)Ft _vG J • (2) 

We neglect the Zeeman splitting, which, in the case of the transverse magnetic field, has a typically much smaller 
effect than the vector-potential term, provided the quasiclassical condition kpl e ^ 1 is satisfied^ 

The Usadel equation is supplemented with the normalization condition for the quasiclassical Green function 

9 2 = l- (3) 

For simplicity, we assume for the moment that the proximity effect is weak (close to the critical temperature of the 
superconductor) and that the boundary conditions at the interface with the superconductor are rigid (which is the 
case for the transparent interface, if the normal region is much more disordered than the superconductor) i 18 ' 19 i 20 We 
will see in Section [V] that these assumptions are not crucial and can be relaxed. 

Then the Green function can be linearized around the normal-metal solution as 

sgnM F \ 
-Ft - B gn(w) J ' [ ' 

and, choosing the gauge A = —yHe x , the linearized Usadel equation (JTJ) takes the forr™^ ,!; 

2 . ^2 M& 



(V, + 2iiryY +vt,- 



hD 

Here we have rescaled both coordinates x and y in the units of the magnetic length 



F(x,y) = 0. (5) 




Zh = \^, (6) 



where 4>q = h/2e is the (superconducting) flux quantum. This equation is supplemented by the boundary conditions 
at the interface with the superconductor and at the open interface, 

F(x = {0,L x },y) = F B e ±4 *, (7) 
V y F(x,y = ±L y /2) = 0. (8) 

The boundary condition |J7J) is the rigid one, with Fg being the bulk value of the anomalous Green function in the 
superconductor (close to the superconducting transition temperature, Fb = A/|w|). The phase difference across the 
junction is thus denoted 2%. Condition (jHJ) expresses the vanishing of the current through the interface with vacuum. 
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FIG. 1: SNS junction in a transverse magnetic field. 

The second anomalous component F^(x,y) is the complex conjugate of F(x,y). The current density can be 
calculated from F(x, y) as¥- 



J = 2meN(0)DT ^ 



F^F - FVFt - l^FFt 



(9) 



where N(0) is the density of states in the normal metal phase (per one spin projection). The symmetry of translation 
along the z direction implies that the current remains in the xy plane. The sum is taken over the Matsubara frequencies 

U). 

III. ANALYTICAL RESULTS FOR A LONG JUNCTION 

We are interested in solving the Usadel equation ([5]) in the middle of the long junction: L x ^> . This will allow us 
to retain only the mode with the slowest decay along the x direction from a spectral decomposition of the solution. To 
simplify our treatment, we assume that the temperature is sufficiently low, compared to the junction length: £t 3> L Xl 
where the thermal length scale is defined as £t = y/hD/T. This assumption allows us to neglect the u> term in the 
Usadel equation ((5]). We will comment on this assumption in Section fVl 

First, notice that in the limit of the long junction (when the junction length L x is much larger than the characteristic 
length of the decay of the anomalous Green function F) we can approximate the solution of as a superposition 21 
of the two Green functions for the semi-infinite SN problem 



F(x, y) w Foo (x, y)e lx + F«, {L x - x, -y)e~ ix 



(10) 



where F ao {x, y) is the solution for the SN problem with the semi-infinite normal layer. It obeys the same equation ([5]) 
with the same boundary condition (jHJ) and with the second boundary condition {7} replaced by F oc (x = 0,y) = Fb 
and F QO (x — * oo,y) = 0. 

It will be convenient to use the Fourier decomposition along the x direction by extending the semi-infinite problem 
to the whole real axis, 



2iiry) 2 + V 2 y 



Foo{x,y) = f(y) 6(x) , 



(11) 



where the right-hand side accounts for the jump in the derivative of the function at x = 0. Taking the Fourier 
transform, we can rewrite this equation in the integral form 



F 00 (x,y)= I ^e lk * V 2 -(fc + 2^) 2 f(y) 



2tt 



(12) 



The function f(y) is fixed self-consistently by the boundary condition Foo^ = 0, y) = 

At positive x, we can close the integration contour in the upper half-plane, and the poles of the integrand are given 
by the zero modes of the operator 



A = V; )- (k + 2th/) 2 



(13) 
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FIG. 2: Effective wave number fc, the lengths are given in units of £h. A purely imaginary A; indicates a monotonic decay. 



(this operator acts on the functions 4>{y) on the interval —L y /2 < y < L y /2 with the boundary conditions i/j'(±L v /2) = 
0). In the long-junction limit, the solution in the middle of the junction is determined by the zero mode with the 
smallest positive imaginary part of k. 

The general solution to the second-order differential equation Alp = can be written in terms of a linear combination 
of two modified Bessel functions,— 



tp = y/k + 2TT V f Cl hi A 

The boundary conditions at y = ±L y /2 fix the ratio C1/C2 and limit the possible values of k to a discrete set. 

In Figure [H we plot the value of k with the smallest positive imaginary part as the width of the junction L y 
increases from zero to infinity. In the limit L y — > 0, the 2ny correction in the operator (|13[) may be neglected, and 
the spectrum is composed of the non-degenerate eigenvalues k = inir / L y (the "leading" eigenvalue with the smallest 
imaginary part in the limit L y — > is thus zero). At a small finite L y , the leading eigenvalue also becomes finite, but 
remains purely imaginary. This follows from the combined symmetry of the complex conjugation and the reflection 
y 1 ► —y, which relates the eigenvalues k and —k*. Since the leading eigenvalue is nondegenerate in the limit L y — > 0, 
by continuity it must remain purely imaginary for sufficiently small iy— 

At larger L yi two imaginary eigenvalues may collide and bifurcate to a pair of complex-conjugate eigenvalues. This 
happens at L y = L c « 0.82 (see Fig. [5]). For L y > L c , we must take into account the contributions of the two 
modes (corresponding to the wave vectors k and —k*) since they decay with the same rate (given by the imaginary 
part of k). In the discussion of the wide-junction limit (Section III.B), we will show that those modes correspond to 
solutions localized close to the two edges of the junction y = ±L y /2 (for L y 3> 1). The critical length L c separates the 
regime where the superconducting anomalous Green function F (x, y) decays along the x direction without oscillations 
(narrow junction, purely imaginary k) and the regime where the decay of the Green function is damped oscillatory 
(wide junction, complex k with both real and imaginary parts). 



(k + 2nyY 



4tt 



C2 if 1/4 



(fc + 2-Ky) 1 



4tt 



(14) 



A. Narrow junction limit 



For L y <C 1 (in the units of we expand the exact solution (I14|) in powers of L y and find the wave number k 
solving the equation for the boundary condition (|8| at y = L y /2. This yields the expansion 

ivr / 4tt 2 4 932tt 4 8 7976tt 6 12 \ 
" V3 y V + "63" y + 218295^ + 13752585 L » + ' ' ') (15) 
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FIG. 3: Superconducting pair correlations normalized to their value at the border of the junction for L y — 0.25£h (dash), 
0.75£h (dot), £ H (dash dot) and 2.5£ H (solid). 



To the lowest order in L y , the solution to the Usadel equation does not depend on y. In this limit, one can simply 
average the y 2 term in the Usadel equation (JSJ) and arrive at a pair-breaking terrn ^ 10 ' 15 ! 23 



with 



™VlF(x) = Qu\+2T)F(x). 



T _ _hD^ _ De 2 H 2 L 2 y 



(16) 



(17) 



4 12ft 

This result obviously reproduces the first term in (|15[) . For wider junctions, the dependence of the anomalous Green 
function F(x, y) along the y direction cannot be neglected anymore, but as long as L y remains smaller than L c , the 
function F(x,y) exhibits a monotonic exponential decay along the x direction. 



B. Wide junction limit 

In the limit of a wide junction L y ;g> 1 (as usual, in the units of £#), the solution is determined by the two complex 
conjugate wave numbers k and — k* . We show below that the asymptotic behavior of k (in the units of ) is 

k = irLy + fc res (18) 

where fc ros is the constant term in the expansion in L~ . 

Indeed, in the wide junction limit, each of the two zero modes (solutions to Atp = 0) is localized near one of the 
two edges of the junction and decays quasiclassically towards the other edge. The solution localized near y = —L y /2 
should therefore have the quasiclassical wave vector in the operator (|13[) vanishing in that region, which immediately 
gives the leading asymptotics k « irL y (the solution localized at the opposite edge has k ps —irL y ). 

To get the subleading term fc ros , we consider one of those zero modes (say, the one localized near y = —L y /2). This 
zero mode decays quasiclassically towards the opposite edge of the junction, and with an exponential precision we 
can replace the boundary condition at y = L y /2 by the decaying condition at infinity, ip(y — > oo) = 0. This selects a 
solution from (fT4")) of the form 



ip = \/k + 2ny K 1 / 4 



{k + 2-Kyf 

47T 



(19) 
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Imposing now the boundary condition ip'(—L y /2) = implies an equation on k res : 

(^res) 



K; 



3/4 



4~ 



0. 



A numerical solution to this equation give; 



,30 



fc rc 



-1.68 + 2.32?. 



(20) 



(21) 



We illustrate our calculation in Fig. [31 where we plot the zero modes below and above the transition. Below the 
transition (for L x < L c ), the solution is nondegenerate and symmetric, while above the transition (L x > L c ) the two 
zero modes are pushed towards the edges of the junction. The characteristic size of the region near the edge where 
the proximity correlations are localized are of the order 1 [from the solution (fl~9|) ]. i.e., in the physical units. 



IV. JOSEPHSON CURRENT 



The transition between the two types of behavior of the anomalous Green function, purely decaying and decaying 
with oscillations, may be observed in the critical current of the SNS junction in a magnetic field. Our result on the 
transition is consistent with the previous numerical work o 9 ' 10 , which indicate that the oscillations appear when the 
width of the normal region L y becomes of the order of the magnetic length We show below that the oscillations 
of the critical current in the SNS system are governed by the same wave vector k as the oscillations of the anomalous 
Green function F oa {x 1 y) in the SN system discussed in Section ITTTT 

In the long-junction limit, the anomalous Green function is given by the expression (|10p . and, using the expression 
([9]), one arrives at the sinusoidal current-phase relation^ 



./tot = lc sin 2x , 

where J to t is the total Josephson current (integrated over the y and z directions). 
In the "pure decay" regime (L y < L c ), the asymptotic behavior of F oa (x,y) is 

F^fay) = F B iP(y)e ikx 



(22) 



(23) 



with a purely imaginary k (see Fig. [2]), and ip(y) proportional to the zero mode of the operator (|13p . This results in 
the exponential decay of the critical current as a function of L x : 



I c = 8neN(0)DTL z £ F B \ - 



Ly/2 



LJ2 



(k + 2ny)il) 2 (y)dy 



,-\k\L x 



(24) 



where L z is the dimension of the junction along the z direction. Note that this expression is real [since tf>(y) = 4>*(—y) 
in this regime] and positive (one checks this numerically). 

In the regime of "decaying oscillations" (L y > L c ), the anomalous Green function F ao {x,y) contains contributions 
from two zero modes, 



Foo(x,y) = F B ^( y y kx +r(-y) 



[here ip(y) ^ ^*{—y) are the two zero modes of the operator 
direction, one finds 



— ik* s 



(25) 



T3")) ]. Integrating the critical current along the y 



I c = 8neN(0)DTL z I ^ F B J Im 



Ly/2 



Ly/2 



(k + 2ny)'4) 2 (y)d y e 



(26) 



so that Im k and Re k describe the rates of decay and oscillations of the critical current as a function of L x , respectively. 
Note that in the case of a wide junction, L y 3> £,h, the localization of the superconducting pair correlations at the 
edge of the junction results in the localization of the superconducting current in the same region. 

We sketch the phase diagram of the junction in Fig. [4] in the coordinates L x and L y . In experiments, however, 
one usually varies the external field for a given junction with fixed dimensions. In this setup, the easiest way to 
observe a transition between the two regimes is to study a junction with L x > L y (theoretically, we assume L x 3> L y , 
but in practice L x may be limited by the thermal length £r and by the smallest measurable critical current). In 
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FIG. 4: Phase diagram of the junction in the magnetic field. The left region corresponds to the pure decay regime while in the 
right one the critical current I c exhibits interference patterns as a function of the field H . 



this case, as the field H increases, one should be able to observe a crossover between the pure-decay regime and the 
decaying-oscillating regime as crosses over L y . At low fields (for 3> L y ), in the pure-decay regime, the field 
dependence of the critical current is 



h— exp 



V3 



(27) 



where I\ is of the order of the critical current in the absence of the field [and we have used the leading term in the 
asymptotics (|15j)]. This expression reproduces the existing result for a SNS junction with a finite depairing rate (|17p 
[in our treatment, the approximation (|10[) implies assuming (j> 4>o)r^ 



h = h- 



V3 <Po 



sinh -£= -f 



(28) 



with I the critical current in the absence of the field. Note that equation (j2"8"| is valid only for linearized Usadel 
equations while we will show in Section V that the domain of validity of the asymptotics (|27[) can be extended to 
non-linear situations. At high fields (for £g <C L y ), the critical current exhibits the decaying-oscillating behavior with 



I c = h-r- ex P 



-2.32-^ 



sm 



— - 1.68-^ 

00 ?H 



^0 



(29) 



HL^Ly is the total 



Here I 2 is of the same order as I\ and Jo (the current in the absence of the external field), < 
flux through the junction, and <po is a phase shift, which we do not compute here. Note that while both expressions 
(j27|) and ((29)) decay exponentially with increasing the field, the expression in the exponent of (|27| is proportional to 
H, while that in the exponent of (|29|) only to \/~H. 

If, however, one considers the current-field dependence for a contact with L x < L y , then one would observe a 
crossover from the Fraunhofer pattern (for 3> L x ^ 10 ^ 15 



In 



sin(7r0/^ o ) 



(30) 



directly to the wide-junction regime (|29p. as the magnetic length becomes shorter than L x 



V. APPLICABILITY OF THE RESULTS 



To simplify our discussion, we have considered in Sections UlrllVI the linearized problem with rigid boundary con- 
ditions. However, our method is based on finding the zero modes of the operator (|13[) which describes the proximity 
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effect in the middle of the junction. Therefore, our results remain valid for more general boundary conditions and 
for the non-linear case, as long as the junction is sufficiently long (so that the Green function F ao (x, y) decays by a 
factor much larger than one by the middle of the junction). In this case, the Usadel equation close to the middle 
of the junction can be linearized anyway, and our treatment of Section [HII can be performed in a similar way, albeit 
with more complicated boundary conditions. Therefore all the conclusions of Sections IIIII and HVl about the different 
interference patterns of the critical current and about the transition value L c rs 0.82^ remain valid. The crucial 
condition for applicability of our method is thus that the junction is much longer than [in the narrow-junction 
limit, we also need to assume that 4> ^ 4>v for applicability of our asymptotic formula (|2"T)) ]. The only role of the 
boundary conditions and of the non-linearity of Usadel equations is the renormalization of the overall coefficients in 
the asymptotics (|27|) and (f2U)) . 

Another approximation used in our calculation is the assumption of low temperature. As we have seen in the 
previous sections, the characteristic scale at which the anomalous Green function F(x,y) varies is of the order £g. 
Therefore, the assumption of low temperature [neglecting uj term in the Usadel equation ((Sj)] implies £t 3> Under 
this low-temperature assumption, the temperature enters only as small corrections to the calculations in the previous 
sections (including corrections to the zero-mode wave vector k). 

We can now qualitatively discuss the high-temperature regime £t < £,h- In this case, the decay of the anomalous 
Green function F(x,y) along the x direction is determined mostly by £r, rather than and, as a result, the 
critical current contains an exponential factor exp(— nL x ) with k v^C-p 1 - However, one can still repeat much of 
the reasoning of section III in the presence of the u term. One then finds that in the limit L x — > oo, even at high 
temperature, one can still distinguish the two regimes of the purely decaying and oscillatory-decaying L x dependence. 
The critical width L c separating the two regimes is slightly decreased as compared to the low-temperature case: a 

i. • . . 2/31/3 

dimensional estimate gives L c oc £jy £t a * ?t *C £h- On the other hand, the same dimensional estimate indicates 

]/3 4/3 

that the field contribution to the decay rate along the x direction is of the order £ T ' £ H , which translates to a 

crossover from the purely Fraunhofer regime to the oscillatory-decaying regime at L x ~ Cj/ 3 ^ 1 ^ 3 , i- e - f° r slightly 
thicker junctions than at low temperatures. 

VI. SUMMARY 

To summarize, we consider a long diffusive SNS junction in an ex ternal magnetic field H. We show that depending 
on the width of the junction relative to the magnetic length £h = \J 4>q/ H two different regimes can be observed. For 
narrow junctions the anomalous Green function F decays monotonically along the junction while for wide junctions 
exponentially damped oscillations are present. We find that the transition between the two regimes occurs for a 
width L c w 0.82£#. Those different behaviors translate in a monotonic decay of the Josephson critical current (|24|) 
as a function of the magnetic flux through a narrow junction and in damped current oscillations for wide junctions. 
Finally, we show that for wide links the current and the superconducting pair correlations are concentrated in a small 
region of size close to the border of the junction. 

The main finding of the present work, in comparison with previous studies of this problem, is the identification of 
the damped-oscillating phase for wide and long junctions. This phase resembles both the Fraunhofer regime (for wide 
and short junctions) and the damped phase (for narrow and long junctions). The period of oscillations is the same 
as in the Fraunhofer interference pattern, while the exponentially decaying factor resembles the damped phase. 

Conceptually, the transition between the two asymptotic regimes for long junctions in our problem is similar to the 
transition between the two regimes in superconductor-ferromagnet-superconductor junctions with domains studied 
in one of our earlier papers^ In both systems, the transition between the purely damped and damped-oscillating 
behavior is related to a bifurcation of the solution to the linearized Usadel equations. 

Experimentally the limit of a long junction is accessible and has been the subject of recent experiments. 15 While the 
decaying regime has been observed, even though without a good quantitative agreement with the theory, the regime 
of decaying oscillations predicted in our present paper has not been reached, because the fields were not sufficiently 
high. In future experimental studies, it may be convenient to use junctions with the aspect ratio L x /L y ~ 1 to access 
this new damped-oscillating regime, in order to be able to use lower fields than for L x /L y 3> 1 junctions. In any 
case, an accurate analysis would be required to distinguish the decaying exponential regime predicted in our paper 
from the distorted Fraunhofer pattern due to possible structural inhomogeneities of the critical current, as discussed 
in Ref. []J. 
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